Guiding-fields for phase-separation: Controlling Liesegang patterns 
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Liesegang patterns emerge from precipitation processes and may be used to build bulk structures 
at submicron lengthscales. Thus they have significant potential for technological applications pro- 
vided adequate methods of control can be devised. Here we describe a simple, physically realizable 
pattern-control based on the notion of driven precipitation meaning that the phase-separation is 
governed by a guiding field such as, for example, a temperature or a pH field. The phase-separation 
is modeled through a non-autonomous Cahn-Hilliard equation whose spinodal is determined by the 
evolving guiding field. Control over the dynamics of the spinodal gives control over the velocity 
of the instability front which separates the stable and unstable regions of the system. Since the 
wavelength of the pattern is largely determined by this velocity, the distance between successive 
precipitation bands becomes controllable. We demonstrate the above ideas by numerical studies 
of a ID system with diffusive guiding field. We find that the results can be accurately described 
by employing a linear stability analysis (pulled-front theory) for determining the velocity - local- 
wavelength relationship. From the perspective of the Liesegang theory, our results indicate that the 
so-called revert patterns may be naturally generated by diffusive guiding fields. 

PACS numbers: 05.70.Ln, 64.60.My, 85.40.Hp, 82.20-w 



I. INTRODUCTION 

Pattern formation is an ubiquitous phenomenon in 
out-of-equilibrium systems, and ordered structures of- 
ten emerge in the wake of a moving reaction front 
There has been recently increasing interest, both exper- 
imental and theoretical, in the study of various types 
of chemically-generated patterns. The main reason is 
that they are expected to provide new bottom-up, self- 
assembling technologies for engineering bulk patterns on 
mesoscopic and microscopic scales (for illustration see 
e -g- [2L y, \% la, la Lfl] from a rapidly-growing bibliography), 
for which the traditional top-down methods (i.e., remov- 
ing material in order to create a structure) are reaching 
the limits of their capabilities. 

Detailed understanding of the mechanisms responsible 
for pattern formation is a key-element in developing tech- 
nological applications since it helps constructing the ap- 
propriate tools for the control of the characteristics of the 
emerging patterns. In this paper we shall focus on design- 
ing a simple method to control the so-called Liesegang 
structures [H, [^]. In particular, this method should be 
useful for a recently-proposed experimental set-up that 
allows to create stamps of such structures 0] . 

Depending on the geometry, Liesegang precipitation 
patterns are bands (in an axially-symmetric configura- 
tion), rings or shells (in a circular, respectively spherical- 
symmetric configuration), clearly separated in the di- 
rection of motion of a chemical reaction front. Several 
generic experimental laws characterize the patterns, see 
e.g. [ToL rill] for reviews. In particular, it is found that the 
positions of the bands usually obey simple laws; e.g. they 
form, with a good approximation, a geometric series with 
increasing distance between consecutive bands. This is 



the so-called regular banding situation, which has been 
recently explained [12| using the phase separation in the 
presence of a moving front as the underlying mechanism. 
Briefly, the reaction front, which moves diffusively, leaves 
behind a constant concentration cq of the reaction prod- 
uct, that we shall conventionally name hereafter C parti- 
cles [H, EH EH • At a coarse-grained level, the dynamics 
of the C particles (that can diffuse, and are also attract- 
ing each other) can be described by a Cahn-Hilliard equa- 
tion 0, E3) EH with a source term corresponding to the 
moving reaction front. Starting with a system free of C's, 
the dynamics of the front brings locally the system across 
the spinodal line, provided that cq is inside the unstable 
region of the phase diagram. A phase separation takes 
then place on a short time scale and a band of precipitate 
is rapidly formed just behind the front. This band acts 
as a sink for the C particles. Then the local concentra- 
tion of C's decreases, bringing the system locally in the 
stable phase again. Thus Liesegang patterns are formed 
since the state of the system at the front is locally and 
quasi-periodically driven into the unstable regime. 

The characteristics of these regular patterns can be 
controlled to some extent through an appropriate choice 
of the concentration of the reagents |19| . of the nature 
of the gel that is fillingthe reaction container (2(|, the 
shape of the container [5| , or through an applied electric 

field [UlEiLHiUl. 

The spinodal decomposition scenario has proved its 
power by describing regular patterns and, furthermore, 
by explaining how those patterns can be influenced by 
the concentration of the outer and inner electrolytes and 
by an external electric field. We will show that it can 
be extended to describe other situations, as well. In- 
deed, there is experimental evidence of Liesegang-type 
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precipitation patterns with decreasing distances between 
successive bands which is termed inverse banding. In 
the borderline case between regular and inverse banding, 
the distances between successive bands are constant, a 
situation called equidistant banding [26| . In our attempts 
of describing the above patterns we were lead to a mech- 
anism which may provide a simple, experimentally real- 
izable control tool of the emerging pattern. 

As described in detail in Sec. [Til our proposal is based 
on a phase separation mechanism in a space- and time- 
dependent guiding field, which could represent, for ex- 
ample, a temperature or a pH field. The pattern for- 
mation is thus modeled through a non- autonomous CH 
equation, whose spinodal line is controlled by the guid- 
ing field. Note that the present design of the guiding 
field is different from the homogeneous (overall) cooling 
that was used in most of the previous studies of non- 
autonomous CH models, see e.g. [13, HH [H, H3| • As we 
shall demonstrate, a simple guiding field is sufficient to 
generate crossover between regular and inverse patterns. 
For example, such a guiding field can be a temperature 
field evolving diffusively due to a temperature difference 
at the boundaries of the system, whose characteristics are 
detailed in Sec. Mil The features of the corresponding 
emerging patterns are analyzed in Sec. |lVl As discussed 
in Sec. [V] our numerical findings can be justified by theo- 
retical arguments relating the velocity of the front of the 
guiding field to the pulled-front velocity resulting from a 
linear stability analysis of the phase-separation process. 
Other, more flexible ways to control the patterns are also 
briefly presented in Sec. lVIl Finally, conclusions and per- 
spectives are discussed in Sec. IVIII 



II. THE MODEL 

Let us consider a tube filled with gel, and an initially 
uniform concentration cq of C particles throughout the 
tube. Assuming axial symmetry along the x-axis of the 
tube, we shall consider that the C-particle concentra- 
tion c(x, t) evolves in time according to the Cahn-Hilliard 
(CH) equation in one dimension. After rescaling the 
space and time variables, this equation can be written 
in the following dimcnsionlcss form 



dc{x, t) 
dt 



dx 2 



ec(x, t) — c 3 (x, t) + 



d 2 c(x,t) 
dx 2 



(1) 



with ^ i ^ L, where L is the dimensionless length 
of the tube. Note that we also performed an appropri- 
ate shift and scaling of the concentration that allows us 
to write the CH equation in a form that is more con- 
venient for the exposition of our problem; namely, this 
form is symmetric with respect to the change in sign of c, 
c — c (see e.g. HI] and footnote [H[ for a more de- 
tailed discussion of this point). The shifted and rescaled 
concentration can take both positive and negative val- 
ues, and the stable configurations are symmetric around 
c = 0. 



The parameter e measures the deviation of the temper- 
ature from the critical temperature T c ; it is negative for 
temperatures above T c (for which no phase-separation is 
possible), while it is positive for temperatures T < T c . 
Below the critical temperature, a uniform concentration 
profile Co inside the spinodal decomposition domain, i.e. 
I Co | ^5 c s = is linearly unstable. A small, lo- 

calized perturbation of the concentration can then trig- 
ger a phase separation throughout the system, through 
the amplification of the unstable modes of wavenumbers 
|fc| < y/e. A large body of work (see e.g. [3l|]) has been 
devoted to the study of the phase separation process in 
the simple case of a uniform parameter e taking the same 
value throughout the system. 

Here we shall concentrate on a different situation, 
namely when e is a field, which evolves according to its 
own dynamics. Moreover, it is possible to control its evo- 
lution. For a simple realization of this control consider 
the following situation. Suppose that at a time t = the 
temperature at one end of the tube is lowered and kept 
at a constant value To < T c thereafter. The other end of 
the tube is supposed to be thermally isolated 33]. The 
temperature profile then evolves in time along the tube 
according to the usual Fourier law of heat conduction, 
and so does the related e(x,t) field, 



de(x, t) d 2 s(x, t) 



dt 



dx 2 



(2) 



where D is the dimensionless thermal diffusion coeffi- 
cient. Through an appropriate scaling of the temper- 
ature, the value of e can be set to —I throughout the 
system at t — 0, while e = +latt^0at the left (x = 0) 
end of the tube; at the right (x = L) end of the tube 
there is no heat flow: 

e{x > 0,t = 0) = — 1 , 
e{x = 0,t) = +1 , 



de , 
d~x {x 



L,t) = 0. 



(3) 



These equations @, © define completely the evolution 
of e(x, t) along the tube, from the onset of the cooling 
procedure till reaching the asymptotic uniform profile e = 
+1 throughout the tube [321 ]. 

For the value e = — 1, the uniform concentration profile 
Co is stable, while for e = +1 it tends to phase-separate 
(i.e., the initial concentration |co| < a/1/3). Therefore, 
with the advancing cooling front, the C-particle concen- 
tration becomes locally unstable with respect to phase 
separation. As a consequence, a pattern made of alter- 
nating low- and high-density phases of C appears simul- 
taneously to the propagation of the cooling front along 
the tube. Its properties and characteristics result thus 
from the CH equation (UJ) coupled to the evolution equa- 
tion ([2]) for e(x,t). Appropriate boundary conditions (i) 
guarantee the conservation of C particles inside the tube 
(more precisely, zero-particle fluxes J c at the edges) , and 
(ii) associated to the initial condition, they also ensure 
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the uniqueness of the solution. The boundary conditions 
we used in our numerical discretized procedure amount, 
in the continuum limit, to 



J c (x = and L, t) = , 
d 3 c 

— (x = andi, t) = 0, 



(4) 



with J c (x,t) — d (ec - 



d 2 c/dx 2 ) jdx. Setting 



d 3 c/dx 3 — means that c at the boundaries relaxes to 
c = ±^/e determined by the boundary value of e. More 
detailed considerations, including other types of bound- 
ary conditions and appropriate discretization schemes, 
are discussed e.g. in Refs. [34|. 

The field e(x,t) related to the diffusive temperature 
profile is thus playing the role of a guiding field. One 
can think of, however, other types of fields e(x,t), other 
boundary and initial conditions for an experimental 
setup. As an example, one can assume that a chemi- 
cal agent is diffusing from one reservoir at the x = 
end of the tube, its concentration changes the local pH 
of the system, and thus may drive the C particles to 
phase separation, etc. Accordingly, we call hereafter e 
the guiding field, and thus shall not restrict ourselves to 
the temperature-like interpretation. 



III. CHARACTERISTICS OF THE DIFFUSIVE 
GUIDING FIELD e(x,t) 

During its time evolution, the guiding field e(x, t) will 
modify locally the position of the spin odal line. At a fixed 
time t, the spinodal density c s — \J e(x, i)/Z will reach 
the value |co| at a given point Xf = Xf(t), therefore initi- 
ating locally a phase separation. The point x f (t) defines 
the position of the instability front, which is thus deter- 
mined by the condition e[x = Xf, t) = 3c 2 ,. Behind the 
front, which propagates to the right, the system becomes 
locally unstable, and phase separates into a precipitation 
pattern of alternate high- and low-density regions of C. 

The diffusion equation (J5]) for e(x,t) with the pre- 
scribed boundary and initial conditions ([3]) can be solved 
through a simple Laplace transform method [35]. One 
obtains for Xf an implicit equation comprising an infi- 
nite sum, 



V 



(-!)« 



f-' 2n + 1 



exp 



(2n- 



1) 2 tt 2 ftD 



X cos 



(2n+ 1)tt 



i- x 4- 



tt(1 + 3c 2 ,) 



(5) 



The resulting trajectory of the instability front Xf(t), as 
well as its velocity Vf(t) = dxf(t)/dt for a particular 
choice of Co are represented in Fig. [TJ Note that when 
the spatial, temporal, and velocity variables are rescaled, 
respectively, by L, L 2 /D, and D/L, as indicated on the 
axis of these plots, then the curves for the trajectory and 
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FIG. 1: Upper panel: Time evolution of the front position for 
Co = —0.05 (solid line). The front moves diffusively for small 
t, Xf(t) ~ \/i (dotted line), while the large-time asymptote 
is given by Eq. (|6]) (dashed line). Lower panel: The velocity 
of the instability front »/ as a function of the front position 
Xf (solid line). The short-time asymptote V/(t) ~ l/v 7 * (dot- 
ted line) and large-time behavior given by Eq. Q (dashed 
line) are also displayed, with P denoting the crossover point. 
The scaling of the spatial, temporal, and velocity variables is 
described in the text. 



velocity of the instability front are universal (for a given 
value of Co). 

As can be seen in Fig. [J the front moves diffusively at 
the beginning, and it accelerates past a crossover point P 
where the acceleration is zero. The large-time asymptote 
for the front position can be obtained by keeping only 



4 



the leading n = term in the sum J5J , 



x f (t)nL 1 



— arccos 

7T 



7T(l+3cg) 



■ exp 



4L 2 



(6) 

This approximate expression is valid provided 
L 2 /2tt 2 D = t mm < t < t max = (4L 2 /n 2 D) 1h{8/[tt(1 + 
3cq)]} where t m i„ is the time when the n = 1 term 
in the sum ([5]) becomes negligible with respect to the 
n = term, while i ma x represents a rough estimate of 
the time it takes the instability front to reach the end 
(x = L) of the tube. The function given by Eq. ^ is 
shown in the upper panel of Fig. [TJ and one can see that 
the asymptote is an excellent approximation past the 
crossover point P. 

The corresponding asymptote for the velocity of the 
front has the remarkable property that, when expressed 
in scaled variables and in terms of the position of the 
front, it becomes independent even of the initial concen- 
tration Cq, 



L 7T 7T / x n 

— Vf = — cot — 1 

D 1 2 2 V L 



(7) 



The above expression is displayed in the lower panel of 
Fig. [T] and one notices again that the approximation is 
very good past the crossover point. 
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FIG. 2: Early stage of pattern formation: snapshots of the 
concentration field c(x) (continuous line), guiding field e(x) 
(dashed line), spinodal lines ±^/e(x)/3 (dashed-dotted line) . 
Dotted lines indicate the local equilibria ±yf e{x). The big dot 
represents the position of the instability front. The parame- 
ters are c = -0.05, L = 1000, D = 4, and t = 4000. The 
inset shows various possible profiles of the guiding field, see 
Sec. IVII namely: the usual diffusive configuration (dashed- 
dotted line); a step-like profile (dashed line); and a rigid 
parabola (continuous line). 



IV. RESULTS 

The coupled non-autonomous CH ([1]) and guiding 
field equations have been solved numerically for dif- 
ferent values of the initial density cq, diffusion constant 
D, and length L of the tube. Figure illustrates the 
early stages of the cooling process, with the profiles of the 
concentration c(x), guiding field s(x), and spinodal line 
±c s (x) = ±y/e(x)/3 at a given time t < t max (before the 
instability front reaches the end of the tube). The con- 
centration field inside the high and low-density emerging 
bands relaxes rather rapidly to the instantaneous, local 
equilibrium values ±y/s(x, t), respectively. 

The pattern initiated by the instability front evolves 
afterwards till reaching a stationary profile, made of al- 
ternate regions of c = ±1 and rather sharp interfaces be- 
tween them, throughout the whole tube. Strictly speak- 
ing, this stationary profile is still evolving through coars- 
ening and band coalescence, as predicted e.g. in [3||. 
However (except eventually for some very closely-spaced 
bands, see below the comments on the plug), its char- 
acteristic evolution time is usually well-beyond any rea- 
sonable experimental time [37| ; from a practical point of 
view one can therefore safely assume its stationarity. 

Three typical stationary patterns of the C-particle con- 
centration field are represented in Fig. O 

Before going into a more detailed analysis, let us enu- 
merate some general qualitative features of the emerging 
patterns: 

(i) The total number of bands increases as D increases, 



for fixed cq and L. 

(ii) For fixed D and L, however, the number of bands de- 
creases with increasing |co| (approaching the spinodal). 

(iii) The first part of the pattern displays regular band- 
ing (i.e. increasing distance between consecutive bands), 
while a second part displays inverse banding. It is impor- 
tant to note that the transition from one type of pattern 
to the other is related to the change in the behavior of 
the velocity of the guiding field, namely from the initial 
diffusive-like motion, to the later-time accelerated one, 
see Fig. Q] 

(iv) In some situations, the pattern contains an initial 
plug, i.e. a rather wide initial region of constant concen- 
tration, see e.g. the second panel of Fig. [3] This effect has 
already been encountered in the usual Liesegang-pattern 
formation (nj El ■ The plug may sometimes result from 
the coalescence, on a time scale of the order 0(^1, of a 
certain number of very closely-spaced bands 37[ ■ A plug 
can also form at the end of the pattern, where the bands 
can be again close enough to each other. Contrary to 
the standard Liesegang pattern whose spatial extension 
is only limited by the length of the tube, in our case the 
length of the patterned region can thus be limited by this 
band-coalescence effect. 

Let us consider now the characteristics of the patterns 
from a more quantitative perspective. Figure [4] shows a 
plot of the band positions x n (which are taken, conven- 
tionally, to be the points where c = 0, with an ascendent 
slope, dc(x n )/dx > 0, and are enumerated in the order 
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FIG. 3: The numerical solution c(x) of the non-autonomous 
Cahn-Hilliard equation ^ in the long-time limit t 2> tmax, 
and for system size L = 4000. Upper panel: Co = —0.05 and 
D = 1; middle panel: co = —0.05 and D = 8; and lower panel: 
cq = -0.2 and D = 1. 



of their appearance, starting from the a; = end of the 



tube) as a function of n; different values of the diffu- 
sion constant D were considered, for fixed Cq = —0.05 
and L = 4000. The presence of a large initial plug may 
have some important effects on the n-dependence of x n 
for small n values. Accordingly, a simple and experimen- 
tally measurable functional expression is only expected 
for sufficiently large values of n, precisely as in the case 
of the usual Liesegang patterns. 
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FIG. 4: The position x n of the n-th band with respect to its 
label n for co = —0.05, L = 4000, and different values of D. 

For the initial, regular-banding part of the pattern, if 
enough bands present, one can fit the positions of the 
bands reasonably well with a geometric series, x n ~ 
exp(np), as for a standard Liesegang pattern fill ]. 
This is obviously related to the initial diffusive-like mo- 
tion of the instability front, that does not differ qualita- 
tively from the motion of the reaction front in the usual 
Liesegang configuration However, a power-law fit- 

ting cannot be excluded either, and further detailed work 
meant to clarify this point is in progress and will be pub- 
lished elsewhere. 

For the second, inverse-banding part of the pattern, the 
positions of the bands for large n-s can be fitted equally 
well as x n ~ Inn or with a power law x n ~ rfi , where 
the exponent (3 « 0.2-0.3 is practically independent of 
D. Since the corresponding distance X n = x n+ i — x n 
between consecutive bands behaves like A„ cx n^" 1 , the 
inequality < 1 ensures precisely the inverse-banding 
character of the pattern. 

Figure [5] displays X n — x n +i — x n as a function of 
n. for the same parameter values as in Fig. |U One no- 
tices clearly the initial region of regular banding (with the 
eventual spurious initial plug) and the inverse banding, 
with the final plug. The tails of these plots for large-rt 
values do not allow to discriminate further between the 
two above-suggested fittings of the relation between x n 
and n in the inverse-banding region. 

Finally, in Fig. [S] we plot the width w n of the n-th high- 
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FIG. 5: The distance A„ = x n +i - x n between two successive FIG. 6: The width w n of the n-th band as a function of the 

bands as a function of n, for c = -0.05, L = 4000, and band posi tion x n , for c = -0.05, L = 4000, and different 

different values of D. One can notice outlier points at the values of D The vert i ca l dashed line indicates the position 

beginning and end of the lines. They correspond to the initial of the crossover po int P (see Sec. [TTTJ) . The jumps at the 

and final plug regions. endpoints of the lines are due to the plugs. 



density band as a function of n. It is remarkable that, 
except for a crossover region between direct and inverse 
banding, one can fit throughout, with a good approxima- 
tion: 

w n nWx n + U. (8) 

For the regular-banding region W > 0, and one can eas- 
ily justify this result simply by using mass-conservation 
arguments for the C-particles, as well as the geometric 
progression of band positions, see [13, EH- However, for 
the inverse banding region W < and the approximate 
nature of this relationhip is related to the fact that the 
band positions are not well-fitted by a geometric series 
in this case. Note that on this figure one can clearly see 
that, as already stated above, the transition from regular 
to inverse banding is marked by the crossover point P of 
the motion of the instability front, see Sec. IIII1 

V. THEORETICAL ARGUMENTS 

Our goal is to devise a simple theoretical approach able 
to explain the characteristics of the patterns observed 
in the numerical simulation of the non-autonomous CH 
equation (fTJ), and to be used further on for predic- 
tive purposes. A basic element of our approach is the 
numerical finding that the characteristics of the pat- 
terns are directly related to the motion of the instability 
front. In particular, the local wavelength of the pattern 
A„ = x n +i — x n (see Figs. [SJB]) is related to the velocity 
Vf =Vf(xf) of the front (see the second panel of Fig. [T]) , 
as discussed below. 

Our approach is based on several assumptions, the va- 
lidity of which is verified a-posteriori by comparison of 



the theoretical findings with the results of the numeri- 
cal simulations. Our first hypothesis is that the guiding 
field moves faster than the diffusing C particles; there- 
fore, the phase-separation does not take place ahead the 
instability front, but only behind it. Note, however, that 
the velocity of the front should not be too high either, 
since otherwise our second hypothesis about the quasi- 
stationarity may not be fulfilled. The meaning of the 
second hypothesis is that although t he local v alue of the 
spinodal concentration c s (x,t) — \J e{x, t)/3 evolves in 
time in the wake of the instability front, this evolution 
can be assumed to be slow enough, so that the local insta- 
bility boundaries associated with the spinodal curve are 
in a quasi-stationary state. We assume therefore that the 
onset of the phase-separation instability is pulled by the 
motion of the guiding front, and consequently we can use 
the standard results of the pulled-front theory 38.f39l.l40j 
to establish the characteristics of the emerging pattern. 

Let us recall here the main results of the standard 
theory. Consider an autonomous CH equation (TT]) with 
e = constant throughout the system, and a uniform un- 
stable concentration cq, \cq\ < c s — \feji. A sharply- 
localized perturbation of this state will then evolve into 
an instability front, with a well-defined velocity, leading 
to phase-separation behind it and to the appearance of 
a pattern of well-defined wavelength. Using linear sta- 
bility analysis arguments, one can easily compute both 
the wavelength A* of the most unstable mode and the 
asymptotic velocity v* of the instability front as a func- 
tion of the distance between the initial concentration and 
the spinodal value. Namely, 

A* = W2(V7 + 2) 

3(V7 + 3) 3 / 2 V ; 
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2(^+2) 
3(V7+l) 1 /2 



,3/2 



(10) 



where a = 3(Cj — Cg) = e- 



"34 



Except for the cases when 



one has band-coalescence (coarsening), this wavelength 
provides the wavelength of the asymptotic emerging pat- 
tern. By eliminating the parameter a between these two 
expressions, one obtains a direct relationship between the 
asymptotic wavelength of the pattern and the asymptotic 
velocity of the instability front, 



A* = 



9.642 

(w*)V3 



(11) 



Note however that the relaxation of the system to this 
asymptotic state goes rather slowly, like (1/t), both 
for the wavelength of the pattern and for the velocity 
of the instability front. Moreover, the transient effects 
tend to increase the wave length of the pattern above 
its asymptotic value A*, see [38Ll39Ll40| for further details. 
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Using the above results of the pulled-front theory, we 
make now the Ansatz that the relationship remains 
valid for our non-autonomous CH equation. More pre- 
cisely, we assume that the local wavelength of the pattern 
is determined by the instantaneous/local velocity of the 
instability front as 



A n 



9.642 



[Vf(xn)} 



1/3 ' 



(12) 



The physical picture underlying the above assumption 
is the following. The instantaneous pulled-front veloc- 
ity v* = Vf(t) dictates, see Eq. (|10[) . an instantanous 
value of the parameter a, let us call it CLfif). This 
means that the local concentration in the vicinity of the 
quasi-stationary instability front adjusts rapidly to the 
value Cf(t) corresponding to the parameter a/, namely 
af(t)=e(x = Xf(t),t)-3c 2 f (t). 

The comparison of the theoretical findings based on 
the above Ansatz with the results of the numerical simu- 
lations is displayed in Fig. [71 where the local wavelength 
of the pattern A ra is plotted versus x n for different values 
of D and L. This figure provides a double-check of the 
Ansatz. Namely, 

(i) If Eq. (|12p is valid, then, since Vf(x n ) is universal 
under appropriate scaling of space, time, and velocities 
(according to Sec. IIII[) . then the plots from the numer- 
ical results should merge when applying the rescaling 
A„ — > \ n (D I X) 1 / 3 and x n — > x n /L. This is, indeed, 
the case, as illustrated by both panels of Fig. [7J 

(ii) All the rescaled plots should fit the theoretical for- 
mula (|12|) shown by solid lines in Fig. [Jj 

The agreement between our simple theoretical predic- 
tions and the results of the simulations is surprisingly 
good. The only exceptions are a few outlier points corre- 
sponding, respectively, to the early band formation (ini- 
tial plug) and to the last bands close to the boundary 
(final plug). There is also a systematic initial and final 
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D=16, L=2000 
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FIG. 7: Comparison of the local wavelengths of the pat- 
tern as determined from numerical simulations (symbols) with 
the theoretical results based on Eq. (|12[) (continuous line). 
The wavelengths X„ and the band positions x n , are rescaled 
as discussed in the text. Upper panel: Co = —0.05 and 
L/D = 1000. Lower panel: c = -0.05 and L/D = 125. 
The vertical dashed lines indicate the position of the crossover 
point P where the instability front has a minimal velocity (see 
Sec. ITTTJ . The outlier endpoints are due to the plugs. 



mismatch, that may be due to the high acceleration of 
the instability front at the very beginning and the very 
end of its motion along the tube (see Fig.[T]), and thus to 
the breaking of the quasi-stationarity hypothesis that lies 
at the basis of our Ansatz. Another origin of discrepancy 
can be the dynamics of the guiding field profile, that may, 
under some circumstances, fail the quasi-stationarity hy- 
pothesis. 

We note that the agreement between numerics and the- 
ory is better for large values of the velocity of the insta- 



8 



bility front, i.e., for smaller values of L/D, as shown in 
the upper panel of Fig. as compared to the lower panel. 
This is probably related to a better, respectively worse 
adequacy of the basic hypothesis of a fast-moving front 
as compared to the diffusion of the C particles. Finally, 
the fact that the numerical wavelengths are systemati- 
cally larger that the theoretically-estimated ones may be 
the combined effect of slow relaxation to the asymptotic 
state and the quasi-stationary nature of our configura- 
tion during the onset of the pattern (see the comments 
above on the effects of transients on the wavelength, in 
the autonomous case). 



VI. PATTERN CONTROL 

We address now the problem of controllability of the 
characteristics of the emerging pattern, ft is obvious from 
the above results that both the qualitative (i.e., regular 
or inverse banding) and quantitative features (like total 
length of the pattern, pattern local wavelength, width of 
bands, etc.) can be controlled in the described configu- 
ration through an appropriate choice of the parameters 
L, Co, and (to a less extent, as more difficult to manip- 
ulate) D. Moreover, these results can be described the- 
oretically in the frame of the pulled-front approximation 
thus providing a method for estimating the parameters 
of the patterns. However, this method of control, al- 
though very simple, is somewhat rigid, since the above- 
mentioned control parameters cannot be changed during 
the process, while, ideally, one requires an easily tuned, 
flexible, external tool of control. One can then think 
about moving the tube with the gel (or maybe a thin 
film of gel) in a prescribed temperature profile, with a 
velocity that can be changed at any moment according 
to the needs. One achieves therefore a guiding field e(x, t) 
that can be externally tuned at any moment and point. 

For example, the simplest configuration one can imag- 
ine is an abrupt, step-like temperature profile that moves 
with velocity Vf, such that e(x,t) = — 1 + 2Q(x — Xf(t)) 

(O^i^ L), where 8() designates the Heaviside step 
function and x/(t) is the instantaneous position of the 
step. If the motion is uniform Vf — const., then one 
obtains equidistant banding. If the motion of the step 
is accelerated or decelerated, then the pattern presents 
inverse-banding, respectively regular banding, with char- 
acteristics that depend on the details of Vf = Vf(t). 

Another simple option is to propagate a smooth, given 
temperature profile along the tube, such that e(x,t) — 
F(x — Xf(t)). Now, the characteristics of the emerging 
pattern do not depend only on the velocity Vf(t) of the 
propagating rigid guiding field profile, but also on the 
shape of this profile. 

In order to illustrate these points, we considered, for 
comparison, the emerging pattern in three situations 
(also illustrated in Fig. (2), namely for: 
(i) The diffusive guiding-field profile, as discussed in the 
previous Sections, for a given set of parameters D, L, and 



Co- Recall that the instability front moves with a velocity 
Vf(t) described in Sec. IIII1 

(ii) A step-like profile of the guiding field that moves with 
the same velocity Vf(t); 

(iii) Finally, a parabolic profile of guiding field, e(x,t) — 
[— 1 + A(xf(t) — x — x ) 2 Q(xf(t) — x — ajo)l • O ne has 
e = —1 for x = Xf + xq and the parameter A is de- 
termined such that for x — Xf one has e = 1 — 3cp. As 
before, this rigid profile moves with the same velocity 
«/(*)• 

The results of the numerical simulation are represented 
in Fig. [51 together with the theoretical result based on 
our Ansatz. One can notice that: 

(i) The pattern can be effectively controlled by the pro- 
posed methods. The effects are qualitatively the same as 
for our usual configuration, but these new methods allow 
for more flexible control. 

(ii) The pattern obtained for the parabolic-like profile 
(with xo = 70) is closer to the pattern obtained for our 
usual configuration, as well as to the theoretical predic- 
tions based on the pulled-front approximation; the pat- 
tern obtained for the step- like profile is much different. 
This convincingly illustrates the importance of the quasi- 
stationarity hypothesis for the pulled-front theory. In- 
deed, this basic ingredient is a good approximation both 
for our usual configuration and for the parabolic profile, 
but it is definitely absent in the case of the step-like pro- 
file, for which the associated abrupt jump in the local 
value of e forbids any possibility of quasi-stationarity. 
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FIG. 8: Local wavelengths of the pattern A n versus band 
positions x n , in appropriate rescaled variables. Symbols: nu- 
merical simulations for three different profiles of the guiding 
field (see the text). Continuous line: theoretical calculations 
based on Eq. (|12|l . The parameters are co = —0.05, D = 16, 
and L = 2000. The vertical dashed line indicates the position, 
along the cylinder's axis, of the crossover point P where the 
instability front has a minimal velocity (see Sec. The 
outliers at the beginning and end of the lines correspond to 
the plugs. 
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VII. CONCLUSION 

We have discussed the problem of how to control pre- 
cipitation patterns by bringing a system described by the 
Cahn-Hilliard equation into an unstable state using a pre- 
scribed guiding held. It was shown that simple, physi- 
cally realizable helds, such as a temperature held gen- 
erated by a temperature jump at the boundary, is suf- 
ficient to generate rather complex precipitation patterns 
even in one dimension. The spacing characteristics of 
the patterns were determined numerically for the case of 
a diffusive guiding held, and we developed a quantitative 
theory for explaining the simulation results. The theory 
is based on relating the velocity of the instability front 
generated by the guiding held to the natural, pulled-front 
velocity of the phase-separation process which, in turn, 
controls the lengthscale of the pattern left in the wake of 
the moving front. 

From a theoretical point of view, our results suggest 
that the inverse-banding phenomena observed in some 
Liesegang experiments may have an explanation in terms 
of a diffusive guiding held. This guiding held is perhaps 
not a temperature held, but may be generated by the 
diffusion of some chemical species which do not take part 



in the reactions and the precipitation but may change e.g. 
the local pH value and thus influences the precipitation 
thresholds. 

As far as the technological applications are concerned, 
it appears that the problem of microfabrication of bulk 
structures by chemical reactions and precipitation 0, d, 
BL H fl is just in the first stages of its development. The 
usefulness of this held will be decided on the possibility 
of creating flexible ways to ensure controllability. Our 
results suggest experimentally feasible solutions for the 
control of a particular precipitation process (formation of 
Liesegang bands). Clearly further studies are necessary 
to develop new methods of control and to sort out the 
question of controllability in more complex cases. 
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